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Effective dynamics in Hamiltonian systems with mixed phase space 
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An adequate characterization of the dynamics of Hamiltonian systems at physically relevant scales 
has been largely lacking. Here we investigate this fundamental problem and we show that the finite- 
scale Hamiltonian dynamics is governed by effective dynamical invariants, which are significantly 
different from the dynamical invariants that describe the asymptotic Hamiltonian dynamics. The 
effective invariants depend both on the scale of resolution and the region of the phase space under 
consideration, and they are naturally interpreted within a framework in which the nonhyperbolic 
dynamics of the Hamiltonian system is modeled as a chain of hyperbolic systems. 
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I. INTRODUCTION 

A comprehensive understanding of Hamiltonian dy- 
namics is a long outstanding problem in nonlinear and 
statistical physics, which has important applications in 
various other areas of physics. Typical Hamiltonian 
systems are nonhyperbolic as they exhibit mixed phase 
space with coexisting regular and chaotic regions. Over 
the past years, a number of ground-breaking works 
HHHHHI1HHH have increasingly elucidated 
the asymptotic behavior of such systems and it is now 
well understood that, because of the stickiness due to 
Kolmogorov-Arnold-Moser (KAM) tori, the chaotic dy- 
namics of typical Hamiltonian systems is fundamentally 
different from that of hyperbolic, fully chaotic systems. 
Here "asymptotic" means in the limit of large time scales 
and small length scales. But in realistic situations, the 
time and length scales are limited. In the case of hyper- 
bolic systems, this is not a constraint because the (statis- 
tical) self-similarity of the underlying invariant sets guar- 
antees the fast convergence of the dynamical invariants 
(entropies, Lyapunov exponents, fractal dimensions, es- 
cape rates, etc) and the asymptotic dynamics turns out 
to be a very good approximation of the dynamics at fi- 
nite scales. In nonhyperbolic systems, however, the self- 
similarity is usually lost because the invariant sets are not 
statistically invariant under magnifications. As a result, 
the finite-scale behavior of a Hamiltonian system may 
be fundamentally different from the asymptotic behavior 
considered previously, which is in turn hard to come by 
either numerically |oL 1 10§ or experimentally [Tl| . 

The aim of this paper is to study the dynamics of 
Hamiltonian systems at finite, physically relevant scales. 
To the best of our knowledge, this problem has not been 
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considered before. Herewith we focus on Hamiltonian 
chaotic scattering, which is one of the most prevalent 
manifestations of chaos in o pen sy stems, with examples 
ranging from fluid dynami cs 1 1 Ot 1 1 1| to solid-state physics 
[l2| to general relativity |l3j|. We show that the finite- 
scale dynamics of a Hamiltonian system is characterized 
by effective dynamical invariants (e.g., effective fractal di- 
mension), which: (1) may be significantly different from 
the corresponding invariants of the asymptotic dynamics; 
(2) depend on the resolution but can be regarded as con- 
stants over many decades in a given region of the phase 
space; and (3) may change drastically from one region 
to another of the same dynamically connected (ergodic) 
component. These features are associated with the slow 
and nonuniform convergence of the invariant measure due 
to the breakdown of self-similarity in nonhyperbolic sys- 
tems. To illustrate the mechanism behind the properties 
of the effective invariants, we introduce a simple deter- 
ministic model which we build on the observation that 
a Hamiltonian system can be represented as a chain of 
hyperbolic systems. 

The paper is organized as follows. We start, in Sec. [HJ 
with the analysis of the invariant measure and the outline 
of the transport structures underlying its convergence. 
Our chain model is introduced and analyzed in Sec. II I II 
The effective fractal dimension is defined in Sec. II VI and 
its properties are verified for a specific system in Sec. [V] 
Conclusions are presented in the last section. 



II. INVARIANT MEASURE 

For concreteness, consider a two-dimensional area pre- 
serving map with a major KAM island surrounded by 
a chaotic region. One such map captures all the main 
properties of a wide class of Hamiltonian systems with 
mixed phase space. When the system is open (scatter- 
ing), almost all particles initialized in the chaotic region 
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eventually escape to infinity. We first study this case 
with a diffusive model for the transversal motion close to 
the main KAM island, obtaining an analytical expression 
for the probability density p(x, t) of particles remaining 
in the scattering region at time t and distance x from 
the island [see APPENDIX]. We find that, in the case of 
chaotic scattering, a singularity develops and the invari- 
ant measure, given by lim^oo p(x, t), accumulates on the 
outermost KAM torus of the KAM island [APPENDIX]. 
Physically, this corresponds to the tendency of nonescap- 
ing particles to concentrate around the regular regions. 
Dynamically, the stickiness due to KAM tori underlies 
two major features of Hamiltonian chaotic scattering, 
namely the algebraic decay of the survival probability 
of particles in the scattering region 0, El- ELS and 
the integer dimension of the chaotic saddle Q, and dis- 
tinguishes this phenomenon from the hyperbolic chaotic 
scattering characterized by exponential decay and nonin- 
teger fractal dimension. However, the convergence of the 
measure is rather slow and highly nonuniform, as shown 
in Fig. n f° r typical parameters, which is in sharp con- 
trast with the fast, uniform convergence observed in hy- 
perbolic systems. Our main results are ultimately related 
to this slow and nonuniform convergence of the invariant 
measure. 

Previous works on transport in Hamiltonian systems 
have used stochastic models, where invariant structures 
around KAM islands are smoothened out and the dy- 
namics is given entirely in terms of a diffusion equation 
P, 0| or a set of transition probabilities (Markov chains 
or trees) 

The stochastic approach is suitable to 
describe transport properties (as above), but cannot be 
used to predict the behavior of dynamical invariants such 
as Lyapunov exponents and fractal dimensions. Here we 
adopt a deterministic approach where we use the Cantori 
surrounding the KAM islands to split the nonhyperbolic 
dynamics of the Hamiltonian system into a chain of hy- 
perbolic dynamical systems. 

Cantori are invariant structures that determine the 
transversal transport close to the KAM islands 0, 0- 
There is a hierarchy of infinitely many Cantori around 
each island. Let C\ denote the area of the scattering 
region outside the outermost Cantorus, C2 denote the 
annular area in between the first and second Cantorus, 
and so on. As j is increased, Cj becomes thinner and 
approaches the corresponding island. For simplicity, we 
consider that there is a single island 0] and that, in each 
iteration, a particle in Cj may either move to the outer 
level Cj-i or the inner level Cj+i or stay in the same level 
3]. Let Aj and A+ denote the transition probabilities 
from level j to j — 1 and j + 1, respectively. A particle in 
C\ may also leave the scattering region, and in this case 
we consider that the particle has escaped. The escaping 
region is denoted by Co- The chaotic saddle is expected 
to have points in Cj for all j > 1. It is natural to assume 
that the transition probabilities Aj and A+ are con- 
stant in time. This means that each individual level can 
be regarded as a hyperbolic scattering system, with its 
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FIG. 1: Snapshots of the probability density p as a function 
of x, for p(x, 0) = 5{x — xo), xo = 1, xi = 2, a = 3, and the 
outermost torus of the KAM island at x = [APPENDIX]. 
The time t is indicated in the figure. 



characteristic exponential decay and nonintegcr chaotic 
saddle dimension. Therefore, a nonhyperbolic scattering 
is in many respects similar to a sequence of hyperbolic 
scatterings. 

III. CHAIN MODEL 

We now introduce a simple deterministic model that 
incorporates the above elements and reproduces essential 
features of the Hamiltonian dynamics. Our model is de- 
picted in Fig. |21 and consists of a semi-infinite chain of 
1-dimensional "/\/-shaped" maps, defined as follows: 

M 3 {x) I - t 3 (x- Aj) + 2, 1/0 <x-Aj < 2/0 
[ 6(ar-l) + l, -l/^<x-l<0, 

where > 3 and < At < 1 - 3/0 (j = 1,2,.. .). 
If x falls in the interval 1/0 < x < 1/0 + A~ , where 
Mj is not defined, the "particle" is considered to have 
crossed a Cantorus to the "outer level" j — 1. This in- 
terval is mapped uniformly to [0,1], and the iteration 
proceeds through Mj_i. Symbolically, this is indicated 
by j 3 ~ 1- Similarly, if x falls into 1 — 1/0 — A^ < 
x < 1 — 1/0, wnere A^ = 1 — 3/0 — A~ ; the parti- 
cle goes to the "inner level", and j — * j + 1. Particles 
that reach I/O < x < 1/O + AT are considered to have 
escaped. The domain of Mj is denoted by Ij and is anal- 
ogous to Cj in a Hamiltonian system, where A~ and A^~ 
represent the transition probabilities. The transition rate 
ratios p = A^ / A~ and v = Aj + i/ Aj are taken in the in- 
terval < p < v < 1 and are set to be independent of j, 
where Aj = A^ + Aj . The parameter p is a measure of 
the fraction of particles in a level j that will move to the 
inner level j + 1 when leaving level j, while v is a measure 
of how much longer it takes for the particles in the inner 
level to escape. The nondependence on j corresponds to 
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FIG. 2: Semi-infinite chain of hyperbolic maps Mj, j — 
1,2,... 



the approximate scaling of the Cantori suggested by the 
renormalization theory |3J. Despite the hyperbolicity of 
each map, the entire chain behaves as a nonhyperbolic 
system. For a uniform initial distribution in 7i, it is not 
difficult to show 16] that the number of particles remain- 
ing in the chain after a long time t decays algebraically 
as Q(t) ~ £~ ln W lnl/ , anc i that the initial conditions of 
never escaping particles form a zero Lebesgue measure 
fractal set with box-counting dimension I. However, the 
finite-scale behavior may deviate considerably from these 
asymptotics, as shown in Fig. [3] 

In Fig. 02a) we show the survival probability Q as a 
function of time. For small \x and i/, the curve is com- 
posed of a discrete sequence of exponentials with scaling 
exponents ln(l — A~), which decrease (in absolute value) 
as we go forward in the sequence. The length of each 
exponential segment is of the order of /i in the decay of 
Q and — ln^ in the variation of hit. This striking be- 
havior is related to the time evolution of the density of 
particles inside the chain. This is shown in Fig. |2Ib), 
where we plot the average position (j) of an ensemble of 
particles initialized in I\ (i.e., j = I). The transitions be- 
tween successive exponentials in the decay of Q [Fig.|3{a)] 
match the transitions from a level j to the next in the 
average position of the remaining particles [Fig. Gib)]. 
In a Hamiltonian system, the increase of (j) in time is 
related to the development of the singular invariant mea- 
sure anticipated in our diffusion analysis [see Fig.Q]. The 
piecewise exponential behavior of Q is smoothened out 
for large fi and v > [i [Figs.EIa) and|3Jb)]. 

In Fig. |3Jc) we show the fractal dimension of the set 
of initial conditions of never escaping particles as com- 
puted from the uncertainty algorithm |17| . which con- 



sists in measuring the scaling of the fraction /(e) of e- 
uncertain points (initial points whose escaping time is 
different from the escaping time of points taken e apart). 
The scaling is statistically well defined over decades and 
the exponent a = Aln/(e)/Alne can be computed ac- 
curately. However, the resulting dimension I — a is not 
only significantly smaller than 1 but also depends crit- 
ically on the region L of the phase space where it is 
computed. The convergence of the dimension is indeed 
so slow that it can only be noticed when observed over 
very many decades of resolution, as shown in Fig. IHId) 
where data of Fig. |^c) is plotted over 35 decades! Ini- 
tially smaller, the dimension measured for L = 1% ap- 
proaches the dimension measured for L — I2 as the scale 
e is reduced beyond 1 0~ 15 (i.e. the corresponding curves 
in Fig. IHId) become parallel). As shown in Fig. Eld), 
this behavior is related to a transition in the average in- 
nermost level {jmax} reached by the particles launched 
from e-uncertain points. As e is further reduced, new 
transitions are expected. The dimension measured in be- 
tween transitions is mainly determined by the dimension 
D = ln3/ln£fc, k = (j m ax), of the corresponding element 
of the chain. For given j and e, the measured dimension 
is larger when L is taken in a denser part of the invariant 
set, such as in the subinterval of I\ first mapped into I 2 
[Fig. El c ); diamonds], because {j max ) is larger in these 
regions. In some regions, however, the measured dimen- 
sion is quite different from the asymptotic value even at 
scales as small as e = 1 0~ 30 . This slow convergence of the 
dimension is due to the slow increase of {jmax), which in 
a Hamiltonian system is related to the slow convergence 
of the invariant measure [Fig. [J. The convergence is 
even slower for smaller /i and larger v. Incidentally, the 
experimental measurements of the fractal dimension are 
usually based on scalings over less than two decades 0] . 
Therefore, at realistic scales the dynamics is clearly not 
governed by the asymptotic dynamical invariants. 

IV. EFFECTIVE DYNAMICAL INVARIANTS 

Our results on the chain model motivate us to intro- 
duce the concept of effective dynamical invariants. As a 
specific example, we consider the effective fractal dimen- 
sion, which, for the intersection of a fractal set S with a 
n-dimensional region L, we define as 

^ ,-r \ dlnf(e') 

where f(e') = N(e')/N (s'), and N(s') and N (e') are 
the number of cubes of e dge length s' needed to cover 
S n L and L, respectively []J| ■ We take L to be a generic 
segment of line [i.e., n = 1 in Eq. Q] intersected by S 
on a fractal set. In the limit e — > 0, we recover the usual 
box-counting dimension D = 1 — liuu^o A In /(e)/ A In e 
of the fractal set SDL, which is known to be 1 for all our 
choices of L. However, for any practical purpose, the pa- 
rameter e is limited and cannot be made arbitrarily small 
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FIG. 3: Chain model for fi = 4.1. (a) Survival probability Q 
and (b) average position (J) as a function of time for fi — 0.01 
and v = 0.02 (full line), y. = 0.01 and v = 0.1 (dashed, bot- 
tom), and = 0.08 and y = 0.1 (dashed, top), (c) Fraction 
/(e) of uncertain points as a function of the scale e for points 
taken from L = Ii (circles), L = I2 (squares), and the subin- 
terval L of I\ first mapped into I2 (diamonds), where [s, = 0.01 
and v — 0.1. Circles in (c) are shifted vertically upward for 
clarity, (d) The same as in (c) for e > 10~ 35 and L = I\ (cir- 
cles), L = I2 (squares), and L = I3 (triangles). Dashed line 
(right-side axis): average maximum j of orbits started from 
e-uncertain points, for L = I\. 



(e.g., it cannot be smaller than the size of the particles, 
the resolution of the experiment, and the length scales 
neglected in modeling the system). At scale e the sys- 
tem behaves as if the fractal dimension were D e ff(L;e) 
(therefore "effective" dimension). In particular, the fi- 
nal state sensitivity of particles launched from L, with 
the initial conditions known within accuracy e*, is de- 
termined by D e ff(L;e*) rather than D: as e is variated 
around e* , the fraction of particles whose final state is 
uncertain scales as e 1 ~ D "ff( L ' e ), which is different from 
the prediction e 1 ^ . This is important in this context 
because, as shown in Fig. [3] (where the effective dimen- 
sion is given by 1 — a), the value of D e ff(L;e) may be 
significantly different from the asymptotic value D = 1 
even for unrealistically small e and may also depend on 
the region of the phase space. Similar considerations ap- 
ply to many other invariants as well. 

We now return to the Hamiltonian case. Consider a 
scattering process in which particles are launched from 
a line L transversal to the stable manifold W s of the 
chaotic saddle. Based on the construction suggested by 
the chain model, it is not difficult to see that W s (~l L 
exhibits a hierarchical structure which is not self-similar 
and is composed of infinitely many nested Cantor sets, 
each of which is associated with the dynamics inside one 
of the regions Cj. As a consequence, the effective dimen- 
sion D e ff(L;e) in Hamiltonian systems is expected to 
behave similarly to the effective dimension in the chain 



model [Figs. Efc) andEJd)]. In particular, D e ff(L;e) is 
expected to display a strong dependence on L and a weak 
dependence on e. 



V. NUMERICAL VERIFICATION 

We test our predictions on the area preserving Henon 
map: f(x,y) = (A — y — x 2 , x), where A is the bifurcation 
parameter. In this system, typical points outside KAM 
islands are eventually mapped to infinity. Because of the 
symmetry / _1 = .90/05, where g(x,y) = (y,x), the 
stable and unstable manifolds of the chaotic saddle are 
obtained from each other by exchanging x and y. For 
A = 0.05, the system displays a period-one and a period- 
four major island, as shown in Fig.^Ja). In the same fig- 
ure we also show the complex invariant structure around 
the islands, the stable manifold of the chaotic saddle, and 
three different choices for the line of starting points: a 
large interval away from the islands (L a ), a small subin- 
terval of this interval where the stable manifold appears 
to be denser (Lb), and an interval closer to the islands 
(L c ). The corresponding effective dimensions are com- 
puted for a wide interval of e. The results are shown 
in Fig. Hb): D eff (L a ;e) = 0.84, D eff (L b ;e) = 0.90, 
and D eff (L c ;e) = 0.97 for 10~ 8 < e < 10~ 5 . These 
results agree with our predictions that the effective frac- 
tal dimension has the following properties: D e ff may 
be significantly different from the asymptotic value 1 of 
the fractal dimension; -D e // depends on the resolution 
e but is nearly constant over decades; D e ff depends on 
the region of the phase space under consideration and, in 
particular, is larger in regions closer to the islands and 
in regions where the stable manifold is denser. Similar 
results are expected for any typical Hamiltonian system 
with mixed phase space. 



VI. CONCLUSIONS 

We have shown that the finite-scale dynamics of Hamil- 
tonian systems, relevant for realistic situations, is gov- 
erned by effective dynamical invariants. The effective 
invariants are not only different from the asymptotic in- 
variants but also from the usual hyperbolic invariants 
because they strongly depend on the region of the phase 
space. Our results are generic and expected to meet 
many practical applications. In particular, our results 
are expected to be relevant for fluid flows, where the ad- 
vection dynamics of tracer particles is often Hamiltonian 
[To). In this context, a slow nonuniform convergence of 
effective invariants is expected not only for time-periodic 
flows, capable of holding KAM tori, but also for a wide 
class of time-irregular incompressible flows with nonslip 
obstacles or aperiodically moving vortices. 
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FIG. 4: (a) KAM islands (blank), stable manifold (gray), and 
the lines of initial conditions (Lf, is a subinterval of L a )- (b) 
Effective dimension for L — L a (circles), L = Lf, (squares), 
and L = L c (triangles). The data in (b) are shifted vertically 
for clarity. 
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APPENDIX 

The diffusion model is: dtP(x,t) = d x [x a d x P(x,t)], 
where P is the probability density of all particles, x > 0, 
and a > 2 [lj. The outermost torus of the KAM is- 
land is at x = 0, where the diffusion rate (proportional 
to x a ) vanishes. In a chaotic scattering process the 
initial distribution of particles is localized apart from 
the confining islands. We take P(x, 0) = S(x — x ), 
x > 0, and consider a particle to escape when it 
reaches x > x%. Under the approximation that for large 
x\ the return of particles can be neglected, we disre- 
gard the boundary condition P(xi,t) — and we take 
the solution to be the corresponding Green function: 
P(x,t) = (a- 2)(xx )-^ 2 yy exp(-y 2 - y 2 )Ip{2yy ), 
where y — (a — 2)~ 1 t~ 1 / 2 x~( a ~ 2 ^ 2 , ya is y at x = xq, 
(3 = (a — I) /(a — 2), and Ip is the modified Bessel 
function, which scales as Ip ~ (2yy ) 13 for small 2yy 
For any fixed x > 0, we can show that the distri- 
bution for large t decreases as P(x,t) ~ where 
(3 = (a — 1)/ {a — 2). On the other hand, as shown in 
Ref. pj, the fraction of particles in the interval x < x\ de- 
cays algebraically as Q(t) = J Q 1 P(x,t)dx ~ t _/3 . Com- 
bining these two results, it follows that the normalized 
probability density p(x,t) = P(x,t)/Q(t) decreases as 
p[x,t) ~ t~ x at each fixed x G (Q,Xi) for large enough t 
and diverges arbitrarily close to x — 0. 
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